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ABSTRACT 

The recent observations of an approximately linear relationship between both Be 
and B and iron in metal poor stars has led to a reassessment of the origin of the light 
elements in the early Galaxy. In addition to standard secondary production of BeB, it is 
necessary to introduce a production mechanism which is independent of the interstellar 
metallicity (primary), and in which freshly synthesized C,0 and He are accelerated 
by supernova shock waves. Primary mechanisms are expected to be dominant at low 
metallicity. At metallicities higher than O/H > —1.75, existing data might indicate that 
secondary production is dominant. In this paper, we focus on the secondary process, 
related to the standard galactic cosmic rays, and we examine the cosmic ray energy 
requirements for both present and past epochs. We find the power input to maintain 
the present-day Galactic cosmic ray flux is about 1.5 x 10 41 erg/s = 5 x 10 50 erg/century; 
this estimate includes energy losses from both the escape of high-energy particle and 
ionization losses from low-energy particles. This implies that, if supernovae are the sites 
of cosmic ray acceleration, the fraction of explosion energy going to accelerated particles 
is about ~ 30%, a value which we obtain consistently both from considering the present 
cosmic ray flux and confinement and from the present 9 Be and 6 Li abundances. 

Using the abundances of 9 Be (and 6 Li) in metal-poor halo stars, we extend the 
analysis to show the effect of the interstellar gas mass on the standard galactic cosmic 
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ray energetic constraints on models of Li, Be, and B evolution. The efficiency of the 
beryllium production per erg may be enhanced in the past by a factor of about 10; 
thus the energetic requirement by itself cannot be used to rule out a secondary origin of 
light elements. Only a clear and undisputable observational determination of the O-Fe 
relation in the halo will discriminate between the two processes. 

Subject headings: cosmic rays - nuclear reactions, nucleosynthesis, abundances 

1. Introduction 

Energetic considerations have historically provided a fundamental constraint on and probe of 
the origin of cosmic rays. For example, energetic constraints have been used by many authors (see, 
e.g., the monograph by Brerezinskii et al. (1990) and references therein) to point out that supernova 
explosions are a plausible cosmic ray power source, while other stars are not. An accounting of 
the energy budget have been particularly useful in the case of Li, Be, and B (LiBeB) production. 
Ryter, Reeves, Gradsztajn, k, Audouze (1970) used this type of argument to rule out in situ LiBeB 
production in young stars by flare particles (the "autogenetic" hypothesis), in favor of the present, 
"galactogenetic," hypothesis of LiBeB production by galactic cosmic rays in the interstellar medium 
(Reeves, Fowler, & Hoyle (1970)). 

Recently, an energetics argument has been invoked to discriminate between two possibilities, 
namely, a primary or a secondary origin of light elements in the halo (Ramaty et al. (1997)). The 
appearance of a linear correlation between BeB and Fe in metal poor stars led to the proposal of 
primary mechanisms (different from standard GCRN) for LiBeB production, in which C, O nuclei 
ejected and accelerated by supernovae fragment on H and He at rest in the interstellar medium 
(e.g., Casse et al. (1995), Ramaty et al. (1998) Vangioni-Flam et al. (1998); Higdon et al. (1998); 
Ramaty k, Lingenfelter (1999)). 

Given the approximately linear correlation between [Be/H] and [Fe/H], and a mean value of 
Be/Fe of about 10 -6 , one can conclude that if 0.1 Mq of iron is produced per supernova, then about 
2 x 10 48 atoms of Be are produced per supernova. Standard secondary models of Be production 
yield an energy requirement of about ~ 5 x 10 3 ergs per Be atom at [Fe/H] = —2.5 or about 10 52 
ergs per supernova going into the production of Be. Of course this number is in excess of the 
total energy available in the supernova explosion. This is the basis of the energetics problem for 
secondary production and will be the main focus of this paper. 

Note that the progenitor of BeB is not Fe but O (and C). Thus the pertinent correlation is 
not between BeB and Fe but between BeB and O. If [O/Fe] is constant then the linearity between 
BeB and Fe translates directly into a linear relation between BeB and O. In this case a primary 
mechanism is clearly required. However, recent observations (Israelian et al. (1998), Boesgaard & 
King (1993), Takeda et al. (2000)), indicate that [O/H] increases with decreasing [Fe/H]. In this 



-3- 



case the linearity between BeB and Fe leads to the correlation: log(Be/H) proportional to about 
1.7 [O/H], and should be explained by a combination of primary and secondary mechanisms (Fields 
and Olive 1999a; Fields et al. 2000, hereafter FOVC). While the new O/Fe results are increasing in 
number they are still controversial. An unambiguous determination of the true O/Fe trend in Pop 
II is crucial, and the results of the present study will only further underscore the need to resolve 
this issue. 

It has been argued (Ramaty et al. (1997), Ramaty, Scully, Lingenfelter, & Kozlovsky (2000), 
Parizot (2000)) that energetic arguments can be used to clarify the primary vs secondary source for 
BeB. These authors have concluded that the secondary process in the early Galaxy meet with strong 
energetic difficulties. However, the energetics computation relies on many relatively uncertain 
factors and deserves a careful analysis (including the sensitivity to the injection energy spectrum 
to test the sensitivity of this criterion). Here, we will consider the energetics of standard Galactic 
cosmic rays in both the early and present Galaxy. 

As we will see, a careful treatment of energetics allows us (1) to clarify and make precise the 
manner in which standard cosmic rays are to be included into chemical evolution, and then (2) to 
examine the viability of proposed cosmic ray origins. We begin by establishing (in §2) the present- 
day cosmic ray power budget, based on the observed properties of Galactic cosmic rays and of the 
Milky Way. By explicitly identifying supernova remnants as the power sources for the cosmic rays 
(in §3), we can express the present day energetics in terms of a cosmic ray acceleration efficiency 
per supernova (in §4). We find that this efficiency is subject to considerable uncertainties, due to 
uncertainties in the parameters describing cosmic rays, the Galaxy, and supernova blast energies. 
Nevertheless, the present efficiency must serve as a fiducial point against which we compare the 
past energetics as determined by Be (§5) and 6 Li (§6). In calculating the energetics of 6 Li and 
Be, we note the importance of the larger gas mass in the early Galaxy relative to today, which 
implies a larger total number of atoms of all kinds in the interstellar medium. This in turn implies 
a larger number of target atoms (for a given metallicity) and thus a larger Be production rate, than 
one would estimate using the present day Galactic gas mass; the upshot is that one finds a higher 
Be production efficiency by a factor of M gaSi i n it/M gaSi today ~ 10 when one takes into account this 
gas-depletion effect (heretofore not considered in previous work). With this factor included, we find 
that the cosmic ray energetics implied by the Be data is in good agreement with present-day values. 
Indeed, we find that the energetics constraint is completely correlated with the basic constraint of 
agreement with Be-Fe data: models which fit the data also do not require more energy efficiency 
in the past. We discuss the implications of our results in §7. 

Thus, the existence of an energetics problem will ultimately rest with new observations of Be at 
low metallicity. A dispersion of points above the secondary production predictions would indicate 
the need for a primary source of Be (perhaps occurring in spatially localized regions) to explain the 
high points. Barring that, the only true discriminant between primary and secondary production 
of Be at this time, will be a final and consistent determination of the O/Fe ratio at low metallicity. 
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2. Standard Cosmic Rays: Formalism and Data 

The energy budget of Galactic cosmic rays is intimately connected with their acceleration, 
propagation, and losses in the Galaxy. Thus, we will follow the sources and sinks of the particles 
and the corresponding energy inputs and losses. The sources of galactic cosmic rays are the events 
which accelerate them to their high energies. In this section, we will simply include the sources 
formally, saving comment on the physical acceleration mechanism (namely supernovae) until the 
next section. The sinks for standard cosmic rays are the losses due to escape from the Galaxy or 
energy deposition (ionization losses) in the ISM. These sources and sinks are described within a 
propagation model. The full, most realistic propagation equation is diffusive, following cosmic rays 
in space as well as time. However, for our purposes we will only need to follow the time history 
of the flux averaged over the Galaxy. Thus, we turn to the simplified propagation scheme of the 
"leaky box" model. Further refinements are in any case unwarranted by the approximations made 
in the chemical evolution portion of the analysis. 



2.1. Cosmic Ray Formalism 

To determine the propagated, interstellar particle spectrum we must first specify a (given) 
source spectrum. The source qE is defined so that the number of new particles (integrated over Att 
steradians) which are injected and accelerated to high energy in time interval dt, volume dV, and 
energy (per nucleon) range (E, E + dE) is dA/jnj = qE dE dt dV; thus qE is the number of particles 
accelerated per unit time in a unit volume and energy band. Similarly, the propagated particles have 
a spectrum Ne defined by their number in a volume dV and energy band dE: dAf pTop = Ne dE dV. 
In other words, Ne is the cosmic ray number density in energy band (E, E + dE). 

In the leaky box model, the source spectrum qE and the propagated spectrum Ne are related 

by 

dN P d 1 

where the b = —dE/dt accounts for ionization energy losses and is defined so that b > 0; for the 
numerical work below we will use the tabulations of Northcliffe & Schilling (1970) and of Janni 
(1982). The effective timescale for inelastic nuclear losses and escape is 1/r = l/r eS c + V r nuc — 
l/r esc . Equation (1) holds for each cosmic ray species; of these, protons and a-particles dominate 
the energy budget. 

We will assume a steady state, dNE/dt = 0, which is valid as long as qE changes slowly 
with respect to the loss timescales; this is an excellent approximation for all but the very earliest 
moments in the Galaxy. In the steady state approximation, eq. (1) is the basic relation used to 
solve for the propagated spectrum given the source spectrum and an escape time. It will also be 
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useful to express Ne in terms of the cosmic ray flux 4>e = Nev, which gives 

where u> = b/pv and the cosmic ray mean path is A = pvT csc ; p is the mean density encountered by 
the cosmic rays. 1 

From these definitions, it follows that the cosmic ray injection power per unit volume is 

/•oo 

{Eq)=A dE Eq E (3) 
J 

where the mass number A accounts for the fact that E is energy per nucleon. The total cosmic ray 
injection power is thus 

W CK ={Eq)V= M (4) 

where V is the volume in which the cosmic rays propagate. M = pV, the gas mass contained in 
V, and since most of the Galaxy's mass today is in the disk, M ~ M g . In fact, Wcr is the sum of 
these terms for each species in the cosmic rays. 



2.2. Cosmic Ray Data and Energy Requirements 

We can now estimate the energy input numerically for the present-day Galaxy, by evaluating 
each of the terms on the right side of eq. (4). Of these, the factors due to cosmic rays, (E q)/p, 
depends on the mean density and on the source spectrum we adopt. The source spectrum is itself 
determined from the observed spectrum, which is fairly well-determined at high energies by the 
local cosmic ray spectrum, but at low (< 100 MeV/n) energies requires extrapolation and thus is 
model-dependent. 

We adopt a source with a momentum (per nucleon) spectrum, q p oc p~ s e~ E ( p ^ Ecut . This has 
the form of a power law in p, characteristic of shock acceleration, with a slope 2 < s < 3. We choose 
a high-energy cutoff E cnt = 10 TeV to account for the maximum acceleration energy of a supernova 
(e.g., Ramaty, Scully, Lingenfelter, & Kozlovsky (2000)). As long as E cut S> 1 GeV, the integration 
of the spectrum and hence the value of Wqr through eqs. (3) and (4) is not very sensitive to the 
choice of E cut . That is, most of the power in the integration comes at lower energies. To convert 



lr rhe mean density can be determined today by comparing A (as measured by the ratio of cosmic ray secondaries 
to primaries, e.g., B/C) and the mean lifetime of long-lived radioactive species in the cosmic rays (such as 10 Be). 
The result (e.g., Simpson & Garcia-Munoz (1988); Lukasiak, Ferrando, McDonald & Webber (1994); Connell (1998)) 
is that p ~ 0.3pg, i.e., cosmic rays encounter a mean density about 1/3 of that of the interstellar medium in the disk 
(rig ~ 1 cm -3 ). This is usually interpreted to mean that the cosmic rays spend a significant part of their lives outside 
of the disk and in the halo, where the gas density is much lower. 



-6- 



this to a kinetic energy spectrum we have q E = q p dp/dE = q p /v and thus 

q E = A v(E)- 1 p(E)~ s e-^/^ut = A (E + m) (E 2 + 2mE)-^l 2 e ~ E l E ^ (5) 

in units where c = 1. This spectrum scales as £ ,_ ( s + 1 )/ 2 in the non-relativistic regime, and as E~ s 
in the relativistic regime. The constants A and s are set by noting that when energy losses can be 
neglected (E 3> 100 MeV/n), eq. (2) gives qE = <PeP/ 'A. This solution, along with measurements 
of 4>e and A, fixes A and s. 

To calculate the integrated spectrum, we must therefore first determine the parameters in eq. 
(5), A and s, from observations of 4>e and A. We first turn to the normalization A. Based on the 
compilation of Mori (1997), we take (f> El = 2.8 cm~ 2 s _1 GcV -1 at the (arbitrary) normalization 
energy E\ = 1 GeV. For the escape pathlength, we adopt the fit of Garcia-Muhoz et al. (1987), 
which has = 9.2 g cm' 2 . This determines the normalization constant in eq. (5) to be A ~ 
E%<f> El /A El . 

The source spectral index s is determined via the measured values of (f> E and A(E); at high 
energies, each of these is found to behave as a power law, and thus q E ~ (f) E /A(E) fixes s. Since 
A(E) is itself determined by measured cosmic ray fluxes (e.g., the B/C ratio), the source spectral 
index is best derived from a consistent analysis which yields both s and A(E). Several such analyses 
have been performed, giving values for s ranging from 2.15 (Buckley, Dwyer, Miiller, Swordy, & 
Tang (1994)) to 2.23 (Engelmann, Ferrando, Soutoul, Gorest & Juliusson (1990)) to 2.3 (Webber 
(1998)) to even 2.6 (Garcia-Muhoz, Simpson, Guzik, Wefel, & Margolis (1987)). Here we wil adopt 
s = 2.2 as a fiducial value, bearing in mind the allowed variation is at least ±0.1. 

Uncertainties are present in both the flux and the escape pathlength. The parameters for the 
flux are sensitive to the extrapolation from the observed, solar system intensity to the interstellar 
intensity. This extrapolation- "demodulation" -becomes important at E < 1 GeV, and is model- 
dependent. Consequently, the combination of measurement errors and demodualtion procedure 
uncertainties leads to variation of about a factor of 5 (see, e.g., Mori (1997)) in the quoted interstellar 
proton flux at 1 GeV, (ft El , and uncertainties in the interstellar spectral index of order ±0.05 units. 
For example, the analysis of Webber (1998) gives (j) El = 5 cm~ 2 s _1 GeV -1 and uses a confinement 
which scales as PR~°' G (with [5 = v/c and R the rigidity) and has A^ = 10.7 g cm -2 ; for s = 2.1, 
this gives an A about 25% larger than our adopted value. 

Both the non-relativistic and relativistic limits show that the input power integral (E q)/p = 
J dE E qE has its maximum contribution in the E mSuX ~ m regime, and taking the dominant 
contribution from the relativistic limit, we find 

(E q)/p * Y {El/m r 2 4^ - SY^ 1 = 2-1 GeV ^ s" 1 (6) 
s — 2 Ae 1 A El 

where the latter expressions use s = 2.2. The factor 



Y = Y, Ay? K ~ 1 + 4(He/p) CR ^ 1.4 (7) 
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includes the additional contribution due to He (with He/H ~ 0.1), with the mass number A; L 
converting from energy per nucleon to energy, and where yp R = ^e/^e 1S the abundance of species 
i in the cosmic rays. We see that the input power scales as <j)E 1 /h-E 1 , as it should, and the the input 
power density is fixed by local measurements of the cosmic ray flux and escape pathlength. We 
can thus analytically estimate (Eq)/n = m p (Eq)/p = 2.5 x 10~ 19 erg yr -1 . This simple analytical 
estimate is to be compared with the results using a full numerical integration of the source spectrum 
for both protons and He, which gives m p {Eq)/p = 4.1 x 10~ 19 erg yr -1 in agreement with our rough 
analytical estimate. 2 If M ~ M g ~ 10 10 M Q , then this gives 

W C r = 4.9 x 10 48 erg yr" 1 = 4.9 x 10 -3 foe yr" 1 (8) 

or 1.5 x 10 41 erg/s. This is similar to other estimates; e.g., Drury, Markiewicz, & Voelk (1989), find 
Wqr = 1.1 x 10 41 erg/s, using a p~ 2 - 2 spectrum but somewhat different gas mass and confinement. 
We note that, for a p~ 2 - 2 spectrum, about 20% of Wqr comes from source energies below 1 GeV; 
this fraction rapidly increases for a steeper spectra, reaching 45% for p~ 2 . Many of these lower- 
energy particles are stopped via ionization losses to the ISM rather than by loss from the Galaxy. 
The contributions to particles such as these are not counted in many estimates of Wqr which only 
include the contributions due to escape losses, but are important and are included in our results. 

The result for Wcr has significant uncertainties tracing back to the input parameters. As noted 
above, the uncertainties due to cosmic ray absolute flux calibrations, spectral index determinations, 
and demodulation corrections, which easily allow for about a factor of 4 variation in A (Mori (1997)). 
As a result, (Eq)/n can vary by a factor of 3 within the current uncertainties. As an illustration 
of the possible variations, and to allow comparison to previous work, we turn to the spectra of 
Ramaty et al. (2000) and Parizot (2000). The Ramaty spectrum has s = 2.5, E cut = 10 GeV/amu, 
and A(E) = 10 g cm~ 2 . To compare with our results, we choose the same flux normalization as our 
fiducial estimate, 4>e x = 2.8 cm -2 s _1 GeV -1 . In this case, we find (Eq)/n = 1.7 x 10~ 19 erg yr _1 . 
This is about 40% of our result. The Parizot spectrum is qE oc E~ 1 e~ E / Ecut , with E cut varying from 
10 MeV to 500 MeV. This is meant to describe acceleration within superbubble shocks. Again, 
for comparison we normalize to 4>e 1 = 2.8 cm" 2 s -1 GeV -1 . For E cut = 500 MeV/amu, we find 
(Eq)/n = 1.2 x 10 -19 erg yr _1 . which is about 1/3 of our result. On the other hand, the E cut = 10 
MeV/amu is a low-energy spectrum which does not resemble the observed GCR flux and should 
not be used in this context. In particular, fixing the normalization to the observed flux at our 
Ei = 1 GeV implies that the flux at, say, E = 10 MeV is higher by a factor of exp(Ei/E cut ) = e 100 , 
with a similarly large increase in (Eq)/n. 

Finally, a significant uncertainty is the present Galactic gas mass. We have used M gas = 
1O 1O M . Smaller estimates exist: Henderson, Jackson, & Kerr (1982) infer a neutral atomic hydro- 
gen content of Mm = 4.8 x 1O 9 M ; including a helium component with mass fraction ~ 0.25 gives 



2 The main difference comes from the differences in the adopted spectrum (eq. 5): and its relativistic approximation 
as AE~^ used in the analytic estimate. When the two are both normalized to <j>E 1 at E\, the power law approximation 
is lower for all E > Ei, accounting for most of the discrepancy. 
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6.4 x 1O 9 M . The total gas mass should include the molecular and ionized components, which each 
are likely to be at least ~ 25% of the neutral atomic component and thus lead to the estimate we 
use. The neutral atomic mass is inferred from 21-cm radio observations, which are themselves quite 
accurate, but in terms of column densities. To arrive at a mass requires assumptions regarding the 
distribution of gas in the Galaxy, which leads to uncertainties of at least a factor of 2. Including 
the gas mass uncertainty with those in the cosmic ray inputs, the estimate of Wqr can easily vary 
by a factor of 5. 

To summarize: we have quantified the power requirement Wqr, and computed the present-day 
value give the flux at Earth. However we have so far made no statement about the physical nature 
of the energy source. We now turn to this issue. 



3. Supernovae as Sources of Cosmic Rays 

We now invoke supernovae as the source - the acceleration engine - of cosmic rays. This 
connection has both observational and theoretical support. Observations of supernova remnants 
confirm that these environments are excellent particle accelerators. Radio synchotron emission has 
long been observed for supernova remnants and establishes that nonthermal, relativistic electrons 
are accelerated in these environments. The detection of X-ray synchotron radiation (Koyama et 
al. (1995)), and TeV 7-rays (Tanimori et al. (1998)) both confirm that the electrons are accelerated 
up to ~ 100 TeV. In additional, recent 7-ray observation provide the first direct evidence of nucleon 
acceleration to E > 300 MeV/n in supernova remnants (Combi, Romero, & Benaglia (1998)). 
These observations are in broad agreement with the predictions of theoretical models of particle 
acceleration in supernova shocks, and thus a coherent picture is emerging in which it seems that 
supernovae are capable of accelerating particles at least to ~ 100 TeV. 

The net Galactic cosmic ray injection power due to supernovae can be written as 

W S n = Vcr E S n A/sn (9) 

where E'sn is the average mechanical energy output of a supernova, A/sn is the Galactic super- 
nova rate, and tjqr = Eqr/Estss is the cosmic ray acceleration efficiency-the fraction a supernova 
explosion's total mechanical energy which goes into cosmic rays. 

The appearance of the supernova rate in eq. (9) will allows us to insert the cosmic ray energy 
budget-and thus the cosmic ray sources and flux-within the chemical evolution formalism in an 
straightforward way. In chemical evolution, the supernova rate is related in a simple way to the 
fundamental input, the global star formation rate ip. Specifically, 

A/sn = dm ^{m)^(t - r m ) ~ ip(t) / dm £(m) = X s ^/{m) S N (10) 

JS7SI JSN 

where £ (m) is the initial mass function and T m is the age of a star of mass m, and "SN" denotes the 
mass range of supernova progenitors, m > 8 — 1OM . The approximation in the last two expressions 
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holds if t S> T mgN ~ 10 Myr, which is an excellent approximation for all but the very earliest 
moments of Galactic history. The supernova mass fraction is ^sn = Jg^ dm m £(m)/ J dm m £(m) 
and the mean supernova progenitor mass is (m)sN = Jsn dm m £(ra)/ Jg N dm £(m). For a Salpeter 
IMF, Xsn = 0.14, and (m)sN = 19M , and thus the supernova and star formation rates are 
connected by Asn = 7.4x 1O _3 ^/M . Thus, if the Galactic supernova rate is about (3±2) x 10 _2 /yr, 
as in a recent estimate (Dragicevich, Blair, Sz Burman (1999)), then one would infer a star formation 
rate ip = 4 M /yr, in general agreement with other estimates (e.g., Timmes, Diehl, & Hartmann 
(1997) and refs. therein; Scalo (1986) and refs. therein). 

Thus, we have the scaling Wsn oc if), with the explicit relation 

W S N = r/CR E S N X S Nlp/{m)sN (11) 

Inserting the galactic supernova rate into (9) or (11), we have 

I^sn = 3 x IO-^cr foe yr- 1 (12) 

where we use .Esn = 10 51 erg = 1 foe, and JVsn = 3 (100 yr) -1 ; these are good fiducial numbers 
but are each subject to uncertainties of at least a factor of 2. 



4. Cosmic Ray Efficiency in the Present Epoch 

Now we can unambiguously link the cosmic ray and chemical evolution formalisms. The 
connection is made by asserting that supernovae are the cosmic ray accelerators. To implement 
this statement we join eqs (4) and (11): 

W C r = I^sn (13) 

This equation has important implications, which we now explore. We first turn to the specific 
numerical implementation of this statement and estimate the present-day cosmic ray acceleration 
efficiency. We then elaborate on how this statement clarifies the formalism used to date in treat- 
ments of cosmic rays in chemical evolution. 



4.1. The Present-day Cosmic Ray Efficiency 

Equations (8), (12), and (13) together imply that the present-day cosmic ray efficiency is 

7?cr = 0.16 (14) 

i.e., 16%, if the supernova mechanical energy is i^sN = 1 foe. We caution that this result is not 
known to high accuracy, due to the uncertainties in the inputs, notably the supernova energy, rate, 
and present gas mass (which also plague previous estimates). If one allows the input parameters 
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to vary within plausible ranges, t/cr can easily span more than order of magnitude, from < 10% 
to 100% or more. (Of course as ?7cr approaches 1, we have the unphysical result that from the 
point of view of eq. (13), cosmic rays require more than the available mechanical energy in the 
supernova.) 

However, it is encouraging that our fiducial result is reasonable and suggestive when one 
recalls that the cosmic ray energy density in the ISM is in rough equipartition with the thermal 
and magnetic energy. In other words, it is well known that the propagated particles reach an energy 
balance with their environment. Here, we see that even in the acceleration process there is a rough 
equipartition of energy between the accelerated particles and the other components. This empirical 
figure is in line with theoretical estimates of particle acceleration efficiencies in supernova remnants. 
For example, Berezhko &: Volk (1997) recently conclude that the acceleration efficiency could reach 
50 per cent. It is also worth noting that this result is directly tied to the local cosmic ray flux. 
Consequently, a similar result would be obtained in all models which fit the observed properties 
of the comic ray flux and propagation. That is, this result is model-independent as an order of 
magnitude estimate, but the precise value depends on the details of the inputs used. 



4.2. Scaling Laws for Cosmic Ray Flux 

Equation (13) implies that 

q = VCR j^r— A SN -i— x — 77 (15) 

where the mean injection energy (E) in j = A J dE Eq^j / dE depends on the shape of the source 
spectrum. Further, since the integrated cosmic ray flux is proportional to the integrated source 
strength, <j) oc q, we have 

oc t (i 6 ) 

Thus, for a fixed volume, a higher supernova rate means a higher star formation rate and more 
cosmic ray acceleration; this in turn means higher flux (for a fixed confinement). To understand the 
appearance of V in eq. (16), imagine dividing the leaky box volume into two identical subregions; 
each of these will have both a star formation rate and a volume that is half of the larger region (recall 
that ip represents the total star formation rate in M0/unit time and is in effect proportional to V). 
Thus eq. (16) is equivalent to the statement that the star formation rate per unit volume must be 
fixed as is the flux of cosmic rays. Finally, if the leaky box volume is fixed-a good approximation 
today, though it may not be for early Galaxy-then we have (f> oc ip, as is usually assumed in the 
literature. 

It is important to note that in the early Galaxy, which gave rise to Pop II stars, conditions 
were likely different from today. These differences can affect not only cosmic ray injection, but also 
cosmic ray propagation. To see this, note that the cosmic ray flux is related to the injection and 
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confinement roughly by <j) ~ Aq/p, which generalizes eq. (16) to 

oc A^/pV (17) 

Thus, the flux strength depends not only on the supernova rate and the volume V of the confinement 
region, but also on the mean escape length A/ p. This ratio (as well as the volume) could have been 
different in the early Galaxy. For example, it is very likely that the early Galaxy as a whole had a 
larger volume (before it settled to a disk), it seems likely that the cosmic ray confinement volume 
was also larger. The cosmic ray escape time would likely have been different as well. The escape 
losses depend on the interplay between the ISM density and the Galactic magnetic field, and as 
it is unclear how these would have varied, and consequently difficult to see whether A/p would 
have been larger or smaller. Clearly it is difficult to ascribe a realistic uncertainty on this quantity. 
Although such a study would be very useful, and would modify the simple scaling of cosmic ray flux 
with the supernova rate, it is clearly beyond the scope of this paper and we will retain the <f> oc ip 
scaling. Another uncertainty concerns the increased gas mass and thus enchanced total number of 
ISM atoms in the halo phase (see below, §5). It could well be that the stars in which LiBeB is 
observed have been formed in the thick galactic disk. 

Indeed this scaling of 4> oc tp has been made by most previous work early Galactic cosmic 
ray nucleosynthesis. Its common use is evidence for its intuitive appeal. Note that here, we have 
derived this scaling. It is demanded by considerations of energetics together with the explicit 
assumption of supernovae as the acceleration engines of Galactic cosmic rays. Thus we see again 
that considerations of energetics help to clarify issues of cosmic ray origin. 



5. Beryllium and Cosmic Ray Efficiency 

So far, we have concentrated on cosmic ray energetics today. Now we turn to cosmic ray 
energetics in the past. Our approach follows the basic argument of Ramaty et al. (1997; 2000), 
which we extend by including aspects of the early Galactic environment that were not heretofore 
emphasized. (Also, our adopted source spectrum is somewhat different, which can be important in 
the comparison of numerical results.) 

The instantaneous cosmic ray input energy per supernova is 

AW C R WcnAt W C r na . 

e = A . r — = — = — (18) 

AAAsn A/" SN At A/sn 

for each epoch t, where Wcr and A/sn are the aggregate injection energy and supernova number 
over the Galactic history. From (18), we can relate e to the cosmic ray efficiency: e = ^cr-^sn- 
To see this relation more clearly, consider models where supernovae power cosmic ray acceleration. 
Then the total injection power is fixed by the scaling Wcr oc qV oc eip, normalized to the present 
value as calculated above. Since A/sn oc tp as well, eq. (18) follows directly. We will assume that e 
is constant in time, and use the Be-0 data to check this assumption. 



-12- 



The observed Be/O ratio is connected to the energetics via WcR/N Be , the ratio of the Galactic 
cosmic-ray power input to the Be production rate N Be . As is well-described in the literature, the 
total Galactic Be mass changes according to the usual chemical evolution expression 

dM Be /dt = m Be Q Be - X Be E ~ m Be AB c ( 19 ) 

Here, Q Bc is the Be source rate, i.e., the total number of Be atoms produced by cosmic rays 
per second. The mass fraction of Be is X Be , and the ejection function E = J dmtp(t — T m )4>(m) 
quantifies the rate at which gas mass is expelled from dying stars. The last approximate equality 
is valid at early times (metallicities less than solar) when Be astration can be neglected. Thus, at 
early times, WcR/N(Be) ~ Wcr/Qbc- This ratio is just the inverse of the Q/W ratio of Ramaty 
et al. (2000), up to the differences in the observed spectra. The N(Be)/WcR ratio depends on the 
CNO abundances in the ISM since AT (Be) oc O/H. 

At early times, therefore, the Be production rate is just given by the source rate 

dN Bc /dt = Q Bc = Y^ Ni fab) (20) 
ij 

where the sum runs over targets % and CR species j, and Aj is the number of Galactic gas atoms 
in the form of i: Ni = mV . It is important to note that the rij refer to the gas density of in the 
Galaxy proper, i.e., in regions whose material will be incorporated into star forming regions and 
ultimately in Pop II stars; this density and is not identical to the mean density n seen by cosmic 
rays on their excursions in and out of the gas-rich parts of the Galaxy. 

In standard cosmic ray nucleosynthesis, the main target is O, and the main projectiles are 
protons. Furthermore, the cosmic ray species are assumed to occur in with the same spectral 
shapes, only different relative abundances, so that (f> l E = yf K cf^ E . Thus we can rewrite 

dN Bc /dt = Q Bc = ^N H %a (21) 
rl 

where $ p = f dE<f>E is the total integrated flux in protons, and 

a = ^^/Oy? R (22) 

is an abundance- and spectrum- weighted cross section, with Ai 6 CNO. 
Thus we have 

dN Bc /dt ~ (0/H)Ahct</> oc (0/B)N H aip (23) 

that is, the product of the number of targets (i.e., Ao = (O/H) Ah), the appropriate spallation cross 
sections, and the cosmic ray flux. The key feature here is the appearance of Ah, the total number 
of H-atoms in the Galactic gas, which encodes the fact that what matters is the absolute number 
of targets and not just their relative abundance. Note that Ah = M gas (H)/m p = XnM gas /m p , and 
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since Ah varies by only a few percent over the history of the Galaxy, then Ah oc M g 
have 

(0/K)N H a<f> (0/R)N R a 



QBe/W 



ZM, 



H'as 



Thus we 



(24) 



(Eq)pV {E)- m] pV/K 

where the second expression uses the mean injection energy defined after eq. (15) and the approxi- 
mation q ~ </>/A. The last expression uses the chemical evolution language of metallicity Z oc O/H 
and gas mass, and makes the usual assumption that A (and thus p) is constant. 

The appearance of Ah (or equivalently, gas mass) in Q/W has not been heretofore emphasized, 
but plays a key role in the energetics, since the gas content (absolute number of atoms) of the Galaxy 
decreases with time, while the metallicity increases. One must include both effects to fully evaluate 
the energetics. Figure 1 illustrates the effect of the M gas factor; without it (dashed curves) one 
simply has N-Q e /W oc Z, and the plot is a straight line. However, with it, the solid curves show that 
there is a peak where the product of decreasing gas mass and increasing metallicity is maximized. 3 
Moreover, at low metallicity, the curve has the same linear shape, but is increased by the factor 
jWgas,init/Af g as,now ~ 10. This large factor plays a key role in the energetics. 

We can thus make a rough analysis of Be energetics in the early Galaxy. This argument was 
proposed and developed by Ramaty et al. (1997; 2000); the detailed treatment here closely follows 
and extends that of Fields & Olive (1999a) and FOVC. The basic idea is that the abundances of Be 
and O at any epoch are a record of the total cosmic ray and supernova energy inputs, respectively, 
for that epoch. Thus an analysis of the evolution of the Be/O ratio can shed light on the evolution 
of the cosmic ray energy budget. 

The cosmic ray input energy per supernova is thus calculable given (1) the "energy per atom" 
ratio, (2) the observed Be/O ratio, and (3) the oxygen output per supernova. Specifically, the 
following expression relates these quantities: 





" ^CR " 




_ A>(Be)~ 




~N(0) 


A"sn 


A>(Be) 


th 


_A>(0)_ 


obs 


_ A/ S n _ 



(25) 



th 



Note that, strictly speaking, the N(Be)/N(0) term is in fact the ratio of Be and O production 
rates in the theory. However, as long as astration effects are unimportant and the Be and O data 
are related by a power law so Be oc O^ 1300 , then 



A>(Be)/A>(0) = ^Bco(Be/0) obs 



(26) 



In this way, one can use the observed Be trends versus O (or Fe) to constrain the cosmic ray energy 
input per supernova over all epochs for which abundances are available. At low-metallicities, the 
Be-0 trend has w Be o - 1-7 ± 0.2. 



3 In the instantaneous recycling approximation, we have M gas 
Ze~ z,>y , which is maximized at Z = y. 



M, 



gas,init 



e- z/y , with y the yield. Thus Q/W 
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Let us apply this analysis to Be-0 data. As seen in Figure 2, the lowest metallicity point 
is [O/H] = —1.74, for which Be/H = 4.4 x 10~ 14 . Using Figure 1, the energy cost is Nbc/W = 
2 x 10~ 3 atom/erg. The data point has Be/O = 2.8 x 10~ 9 and using an oxygen yield per supernova 
of M(0)/jvsN = 2M Q corresponding to 1.5 xlO 56 atoms of oxygen, we find that we are required to 
produce 4.2 x 10 47 atoms of Be per SN. We thus find that the cosmic ray energy input per supernova 
at this metallicity is e = 0.36 foe/SN, corresponding to an efficiency of t/cr = 0.36 for Esn = 1 
foe. Thus we see that on the basis of Be production alone, one requires a total energy available per 
SN that is at about the level required to produce the observed flux of cosmic rays. In other words, 
we are able to produce the desired abundance of Be at low metallicity with an efficiency which is 
comparable that needed to accelerate cosmic rays. 

The importance of the Ah factor in the energy per atom calculation is now clear, as it lowers 
the cosmic ray energy budget to within values that lie within the wide range allowed by present-day 
data and theory. Without this factor, the energy requirements inferred from Be would increase by 
about an order of magnitude. This would put the needed energy budget beyond the traditional 
estimates. It is worth noting, however, that even this case is not out of reach, if supernova blast 
energies are more typically 3 foe, and the cosmic ray efficiencies approach 50-60%, notwithstanding 
the other uncertainties discussed in §2.2. 

Let us compare the above result with the previous arguments made by Ramaty et al. (2000) 
and Parizot (2000). As mentioned in the introduction, Ramaty et al. base their argument on the 
evolution of Be with respect to Fe. At the low oxygen abundance of [O/H] = -1.7, we can assume a 
value of [Fe/H] = —2.5. Recall, that the conclusion that secondary Be production was energetically 
excluded was based on the Fe yield, the constancy of Be/Fe, and the determination of the number 
of Be atoms per erg produced per supernova. At this metallicity, the Be/O ratio of 2.8 xl0~ 9 
corresponds to a ratio of Be/Fe of 3.2 xl0~ 7 and is significantly below the constant mean value of 
10~ 6 . This is a result of the fact that the data show that the Be/Fe is not a constant with respect 
to Fe/H, but decreases with decreasing metallicity (though slowly) and that there is considerable 
dispersion in the data. Had we used the estimate of energy requirement for secondary production 
of Ramaty et al. (2000) (4 xl0~ 4 atoms/erg at this metallicity), we would have concluded that 
(1.7 x 2.8 x 10~ 9 )(1.5 x 10 56 )/4 x 10~ 4 = 1.7 x 10 51 ergs/SN is required. This may or may not 
be too much depending on the actual energy available. We have typically assumed that the total 
energy available in the explosion is 10 51 erg, but this is uncertain by a factor of at least 3, and if the 
available energy is actually 3 times higher, there is no energetics problem for the lowest metallicity 
points (and lowest Be) points observed. It is however very likely that primary sources are required 
at these low metallicities in order to explain the scatter in the data. Note that as argued above, the 
energy required per Be atom is significantly reduced when density effects are taken into account. 
In that case, energetics is unable to provide a strong indication as to the primary vs secondary 
nature of Be production. 

Parizot (2000) also provided an energetics argument against secondary production. We find 
several differences in the simple estimate which is used to preclude secondary production. First, 
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at low metallicity, Parizot assumes a Be/O ratio of 4 x 10 -9 , where it is in fact 2.8 x 10~ 9 (this 
stems in part from his incorrect value of the solar oxygen abundance) . He further assumes that the 
maximum energy available for Be production is 10 50 erg/SN corresponding to an efficiency of 10%. 
This could be 3-5 times higher and as noted above, the total SN energy could be 3 times higher. 
Finally his choice of a spectrum is far more pessimistic than even that used in Ramaty etal and 
assumes a high energy cut off of 0.5 GeV. Thus from the argument given in Parizot (2000), one 
can not place a constraint on the secondary nature of Be. 

For late times and high metallicities, the approximations of no Be astration and a constant Be- 
O slope break down. A more refined analysis is needed. Given the data, what is actually observable 
is the Be/O ratio, which is in general not set by the ratio of mass change rates 7V(Be)/iV(0), but 
rather Be/O. Since these observables connect not to rates but to integral quantities, they do not 
allow for direct measure of the instantaneous cosmic ray energy input Wqr- Instead, the data can 
determine the time-averaged value 



W C R(t) = Sq dt' W C r = Sj dt' A~sn e 
AM*) So dt ' -^sn Jo dt' Msn 



(27) 



Within the model, we assume (and then test) a constant energy budget e = eo = const, and thus 
(e) = eo as well. However, this quantity can also be calculated in a way that makes explicit use of 
the observed Be abundances. 



We can calculate the (time-averaged) energetics as follows: 
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The cosmic ray input energy per supernova is thus calculable given (1) this "energy per atom" 
ratio, (2) the observed Be/O ratio, and (3) the oxygen output per supernova. In this way, one can 
use the observed Be trends versus O (or Fe) to constrain the cosmic ray energy input per supernova 
over all epochs for which abundances are available. To do this, we rewrite eq. (28), noting that 
in the code all rates in eq. (28) are computed explicitly. Thus we are free to rearrange terms and 
write 
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(29) 
(30) 



where e t h is the theoretical efficiency (and must be less than 30 -50%). Thus the true needed 
efficiency e will exceed the theoretical efficiency only if the observed Be/O ratios are greater than 
the theoretically predicted ones. That is if there is significant (upwards) scatter in the Be vs O 
data, or there is a true excess in Be/H at low metallicity, will we find that energetics preclude the 
secondary production of Be as was argued in FOVC. 
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Thus the energetics entirely reflects the agreement of the abundance trends. If the data are 
fit within the errors, then the right side of eq. (30) is constant in metallicity, consistent with the 
assumption of constant Wcr/JVsn- The is the case for primary models of, e.g., Vangioni-Flam et 
al. (1998) or Ramaty et al. (2000). This agreement with the data also holds case for standard GCR 
models with changing O/Fe, at least down to metallicities [O/H] ~ —1.5. At lower metallicities, 
the data is not conclusive but seems to suggest the need for a primary component (FOVC). If this 
is borne out by improved data, then the need for a primary component of Be will also be reflected 
in a energetics problem for the standard GCR model at these metallicities. Careful observations at 
these metallicities are thus crucial and encouraged. In any case, the standard GCR model certainly 
provides an acceptable fit to the Be-0 data above about [O/H] ~ —1.5 (or [Fe/H] ~ —2.2), and 
thus suffers no energetics problem in this regime. 

To illustrate this conclusion, we can calculate the energetics for all of the data points illustrated 
in Figure 2. We implement eq. (30) for each data point, using our model to calculate [e/(Be/0)] t h, 
and the data for {Be/0) ^ s . The result appears in Figure 3, both with and without normalizing 
the Be solar abundance to its observed value. We see that in the standard case, there is scatter in 
the data but the mean remains close to the solar value, which is around 0.08 foe/SN. 

Combining these factors, the solar beryllium and oxygen abundances give a present-day cosmic 
ray energy budget per supernova which spans a range of 

e = 0.5 ± 0.4 foe (per SN) (31) 

(using the weighted mean and the 2 — a weighted sample variance) or an efficiency of t/cr = (50 ± 
40)% for a supernova blast wave energy of 1 foe. This number is consistent, within uncertainties, 
with our estimate of t/cr = 16% using the observed cosmic ray flux (eq. 14). This concordance 
means that the present-day Be abundance is consistent with the energy that has gone into the 
cosmic rays up to the present epoch. Thus, standard Galactic cosmic rays are able to account 
for both the present-day Be abundance as well as the associated energy, within the described 
uncertainties. 

An important but somewhat technical detail concerns the normalization of the theory curves. 
The approach we have taken in this and previous work is to normalize the curves by demanding 
that the Be abundance at solar metallicity is equal to the solar Be abundance. Physically, this 
amounts to determining the scaling between the global average Galactic flux at the present epoch 
and the present flux at the solar system. The normalization factor is about 0.5, i.e., the solar flux 
is slightly high. This factor has been included in Figure 1. If we did not include this factor, we 
would both raise the entire Be-0 curve and increase the energy requirement by the same factor of 
1/0.5 = 2. The overall mean is thus higher in the non-normalized case, at 1.0 ± 0.8 foe/SN. This 
is just at the fiducial (but poorly determined) supernova mechanical energy budget, but of course 
the lack of normalization also leads to a poorer fit to the Be-0 data. Thus we see an illustration 
of the connection between energetics and the Be-0 fit: a poorer agreement with Be-0 data goes 
hand in hand with poorer agreement with energetics requirements. 
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Finally, we recall that the energy budget we have determined depends sensitively on the cosmic 
ray spectrum adopted. We have used a p~ 2 ' 2 source, which we consider to be a fair representation 
of the experimental determinations; however, the present cosmic ray data allow for other spectral 
forms and different energetic requirements. For example, used by Ramaty et al. (2000) is more 
energetically efficient by a factor of 2.5 (i.e., (Eq) is 40% of the value using our adopted spectrum), 
as noted in §2.2. Thus, using that spectrum, the energy requirements would drop by a factor of 
2.5, both for the present cosmic rays and for the past as derived form Be. Thus, while the ratios 
of past/present values would be similar, but the absolute energy requirements would be reduced. 
Using the Parizot spectrum (with 500 MeV cutoff) would reduce this requirement yet further. 

To summarize, for models with supernova-powered cosmic rays, the fit of the energetics is 
completely determined by, and follows from, the fit of the observed Be-0 data. Models which fit 
the data (and have ip cx 7Vsn) do not have an energetics problem in the past if they do not have 
one now. Conversely, models which do not fit the data (in particular, those which underpredict the 
Be abundance) can also be expected to have energetic difficulties over the Be-0 regimes which are 
poorly fit. 



6. Lithium-6 Constraints on Energetics 

We can repeat this analysis for other cosmic-ray produced elements. For 7 Li and n B, the 
results are complicated by the fact that these elements are produced in other sources. However, 
6 Li, like Be, is a pure cosmic-ray element, and thus the energetics evaluation should proceed as 
it does for Be. We first consider the present-day energetics (the analog of eq. 31). For 6 Li, we 
have Qq^/Wcr = 0.05 atom/erg using the solar value from Figure 1, and ( 6 Li/O) = 1.8 x 10~ 7 . 
Finally, we again take an average of 2M Q of oxygen ejected per supernova, or 1.5 x 10 56 atoms. 
Together, this gives 

e = 0.54 foe (per SN) (32) 

This is comparable to, albeit somewhat higher, than the other estimates we have made. Also, as 
with the Be analysis, using a source spectrum such as that of Ramaty et al. (2000) would lower 
this requirement by a factor of 2.5. 

Regarding the energetics of 6 Li in the past, the analysis of the previous section carries over: 
the efficiency remains constant if the theory and data agree. The data for 6 Li in Pop II stars 
are difficult to obtain, and thus this is a less definitive test. However, as shown by Fields & Olive 
(1999b), the predictions of standard GCR agree with the presently available data. Thus we see that 
the 6 Li constraints on standard GCR are quantitatively similar to those of Be, and demonstrate 
the basic consistency of the analysis. 

Note that 6 Li can be made via the a + a fusion reaction. This production mode is always 
primary since most helium in the ISM is primordial and hence predates the Galaxy. 6 Li is thus 
primary in both the standard GCR scenario of particle acceleration from the ISM, as well as in 
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the case GCR origin is dominated by superbubbles (Ramaty et al. (2000)). An approach which 
synthesizes aspects of these two outlooks adds to the standard GCR scenario a component of 
energetic particles (EPs) which could also originate from superbubbles (for a review see Vangioni- 
Flam, Casse &; Audouze (2000)). The EP component is highly enriched in a particles since its 
origin is in the massive star ejecta which fill the superbubbles. Due to this enhancement, the EP 
component can explain the 6 Li data with an even lower energy efficiency than that required above 
for GCR particles. Further 6 Li data would help distinguish between the EP and GCR scenarios for 
the a-component of early Galactic cosmic rays. 

7. Conclusions 

In this paper, we have reviewed the past and present energy budget for standard Galactic 
cosmic rays We determine present-day cosmic ray power requirements from the observed cosmic 
ray flux and confinement. Comparing this to the fiducial 10 51 erg of kinetic energy in a supernova 
blast, we find the efficiency of cosmic ray acceleration is about 30%. This value cannot be calculated 
to a precision better than a factor of about 3, given the uncertainties in the cosmic ray and Galactic 
input parameters. Nevertheless, we find that this value is consistent with estimates derived from 
present (solar) 9 Be and 6 Li values. 

We have extended our analysis to the study the cosmic ray acceleration efficiency in the early 
Galaxy, as encoded in the Be abundances of metal-poor stars. We find that energy considerations 
are indeed powerful, and allow us to derive the usual <f> oc ip scaling. We also find that in general, 
the requirement that models satisfy energetic requirements is equivalent to the requirement that the 
models produce LiBeB evolution which fits the abundances data, if SN are cosmic ray acceleration 
engines. 

Furthermore, constructing a careful and consistent model for the energetics of cosmic ray 
acceleration and LiBeB production has an additional benefit: it brings to the surface key underlying 
physical issues and assumptions regarding early Galactic LiBeB production. Specifically, we have 
emphasized that the ISM gas content, and thus the total number of ISM atoms, was higher in 
the early Galaxy. Thus, for a fixed metallicity, the total Galactic Be production is enhanced due 
to the increased total number of targets. This effect has been neglected up to now but can play 
a crucial role in the energetics: due to this effect, the number of Be atoms produced per cosmic 
ray erg should be multiplied by about a factor 10 in the early Galaxy compared to the estimates 
lacking this effect. Consequently, the energy test is less stringent than previously thought. Note 
that one could invoke still more complications in the form of early Galactic evolution of cosmic ray 
confinement properties, or a lack of equilibrium between cosmic ray sources and sinks; further study 
of these issues is best done in a consistent model of cosmic ray acceleration and propagation in a 
Galaxy which evolves in time and space. In the absence of such a model, the current "leaky box" 
picture that has been adopted for Be studies, with constant cosmic ray confinement parameters, 
leads to the factor of ~ 10 increase in Be production efficiency we have described. With this effect, 
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the standard GCR model is able both to fit the Be-Fe data, and to maintain a reasonable (and 
roughly constant) acceleration efficiency over the history of the Galaxy at least down to [Fe/H] 
~ -2.5. 

Thus we find that energetic arguments are extremely useful tools in analyzing LiBeB production 
in the early Galaxy, and thereby to examine questions of cosmic ray origin. However, contrary to 
some analyses, we find that the standard GCR scenario does not suffer energetics problems, and 
thus that energetics alone cannot rule this model out. Therefore, the key of the origin of early 
Galactic cosmic rays and LiBeB (i.e., a primary or secondary process?) remains in the hands of 
the observers. Specifically, it is crucial to reconcile the different measurements of oxygen (and thus 
BeB-0 and O-Fe) in halo stars. The determination of reliable oxygen abundances is to be ranked 
as a first priority of the abundance observation programs. 

Indeed, looking towards the future, we can hope to turn the problem around. With reliable 
and accurate oxygen data, we can establish the origin of LiBeB and cosmic rays in the early Galaxy. 
With this knowledge in hand, we can then use the detailed LiBeB abundance patterns and scatter 
to learn more about the physical conditions of gas and accelerated particles in the early Galaxy. 

We thank Tom Jones and Vasiliki Pavlidou for fruitful discussions. This work was supported 
in part by DoE grant DE-FG02-94ER-40823 at the University of Minnesota. 

REFERENCES 

Berezhko, E. G. and Volk, H. J. 1997, Astroparticle Physics, 7, 183 
Boesgaard A. & King J.R. 1993, AJ, 106, 2309 

Brerezinskii, V.S., Bulanov, S.V., Dogiel, V.A., Ginzburg, V.L., Sz Ptuskin, V.S. 1990, Astrophysics 
of Cosmic Rays, (Amsterdam: North-Holland) 

Buckley, J., Dwyer, J., Mueller, D., Swordy, S. and Tang, K. K. 1994, ApJ, 429, 736 

Casse, M., Lehoucq, R., k Vangioni-Flam, E. 1995, Nature, 373, 38 

Connell, J. J. 1998, ApJ, 501, L59 

Combi, J.A., Romero, G.E., & Benaglia, P. 1998, A&A, 333, L91 
Dragicevich, P. M., Blair, D. G., & Burman, R. R. 1999, MNRAS, 302, 693 
Drury, L. O., Markiewicz, W. J. and Voelk, H. J. 1989, A&A, 225, 179 

Engelmann, J. J., Ferrando, P., Soutoul, A., Goret, P. and Juliusson, E. 1990, A&A, 233, 96 
Fields, B.D., & Olive, K.A. 1999a, ApJ, 516, 797 



-20- 



Fields, B.D., & Olive, K.A. 1999b, New Astron., 4, 255 

Fields, B.D., & Olive, K.A., Vangioni-Flam, E. & Casse, M. 2000, ApJ, 540, 930 (FOVC) 

Garcia-Munoz, M., Simpson, J.A., Guzik, T.G., Wefel, J.P., & Margolis, S.H. 1987, ApJS, 64, 269 

Henderson, A.P., Jackson, P.D., & Kerr, F.J. 1982, ApJ, 263, 116 

Higdon, B., Lingenfelter, R.E., & Ramaty, R. 1998, ApJ, 509, L33 

Israelian, G., Garcia-Lopez, R.J., & Rebolo, R. 1998, ApJ, 507, 805 

Janni, J.F. 1982, Atomic Data Nucl. Data Tables, 27, 34 

Koyama, K., et al. 1995, Nature, 378, 255 

Lukasiak, A., Ferrando, P., McDonald, F. B. and Webber, W. R. 1994, ApJ, 423, 426 

Northcliffe, L.C., & Schilling, R.F., Nucl. Data Tables, A7, 233 

Mori, M. 1997, ApJ, 4478, 225 

Parizot, E. 2000, A&A, in press (astro-ph/0006099) 

Ramaty, R., Kozlovsky, B., Lingenfelter, R.E., & Reeves, H. 1997, ApJ, 488, 730 

Ramaty, R., Kozlovsky, b., & Lingenfelter, R. E. 1998, Physics Today, 51, no. 4, 30 

Ramaty, R., & Lingenfelter, R.E. 1999, in LiBeB Cosmic Rays and Related X- and Gamma- Rays, 
eds. Ramaty et al. , ASP, Vol. 171, p. 104 

Ramaty, R., Scully, S.T., Lingenfelter, R.E., & Kozlovsky, B. 2000, ApJ, 534, 747 

Reeves, H., Fowler, W.A., & Hoyle, F. 1970, Nature, 

Ryter, C, Reeves, H., Gradsztajn, E., & Audouze, J. 1970, A&A, 8, 389 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1 

Simpson, J. A. & Garcia-Munoz, M. 1988, Space Sci. Rev., 46, 205 

Takeda, Y., Takada-Hidai, M., Sato, S., Sargent, W.L.W., Lu, L., Barlow, T.A., Jugaku, J. 2000, 
astro-ph/0007007 

Tanimori, T., et al. 1998, ApJ, 497, L25 

Timmes, F. X., Diehl, R. & Hartmann, D. H. 1997, ApJ, 479, 760 
Vangioni-Flam, E., Casse, M., & Audouze, J. 2000, Phys. Reports, 333-334, 365 



- 21 - 



Vangioni-Flam, E., Casse, M., Fields, B.D., & Olive, K.A. 1996, ApJ, 468, 199 
Vangioni-Flam, E., Ramaty, R., Olive, K.A., k Casse, M. 1998, AA 337, 714 
Webber, W.R., 1998 ApJ, 506, 329 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



-22- 



FIGURES 
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1. A plot of Qeb/Wcr for our model. Dotted lines: constant iVn; full lines: evolving N^. 
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2. Be and O data, derived from a consistent set of stellar atmospheres using the "Balmer" 
method described in Fields et al. (2000). The curve is for a model in which Be is produced 
by standard Galactic cosmic rays. 
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3. Energetic requirements inferred point-by-point from the Be-0 data and the model shown in 
Figure (2). The upper panel is for the standard model which normalizes the Be point to solar; 
the lower panel does not use this scaling. 



